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Abstract 

Energy efficiency and its optimization constitute critical tasks in the design of low-power wireless 
networks. The present work is devoted to the error rate analysis and energy-efficiency optimization of 
regenerative cooperative networks in the presence of multipath fading under spatial correlation. To this 
end, exact and asymptotic analytic expressions are firstly derived for the symbol-error-rate of M —ary 
quadrature amplitude and M—ary phase shift keying modulations assuming a dual-hop decode-and- 
forward relay system, spatially correlated Nakagami —m multipath fading and maximum ratio combining. 

The derived expressions are subsequently employed in quantifying the energy consumption of the con¬ 
sidered system, incorporating both transmit energy and the energy consumed by the transceiver circuits, 
as well as in deriving the optimal power allocation formulation for minimizing energy consumption 
under certain quality-of-service requirements. A relatively harsh path-loss model, that also accounts 
for realistic device-to-device communications, is adopted in numerical evaluations and various useful 
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insights are provided for the design of future low-energy wireless networks deployments. Indicatively, 
it is shown that depending on the degree of spatial correlation, severity of fading, transmission distance, 
relay location and power allocation strategy, target performance can be achieved with large overall 
energy reduction compared to direct transmission reference. 

Index Terms 

Energy efficiency, error rate, maximum ratio combining, multipath fading, optimization, power 
allocation, quality-of-service, regenerative relaying, spatial correlation, asymptotic analysis. 

I. Introduction 

Emerging communication systems are expected to provide high-speed data transmission, effi¬ 
cient wireless access, high quality of service (QoS) and reliable network coverage with reduced 
processing time and energy as well as widespread use of smart phones and other intelligent 
mobile devices. However, the currently witnessed scarcity of the two core fundamental resources, 
power and bandwidth, constitutes a significant challenge to satisfy these demands while it is 
known that wireless channel impairments such as multipath fading, shadowing and interference 
degrade information signals during wireless propagation. Furthermore, most energy constrained 
devices, such as terminals of mobile cellular, ad-hoc and wireless sensor networks, are typically 
powered by small batteries where replacement is rather difficult and costly [JT), El- Therefore, 
finding a robust strategy for energy efficient transmission and minimized energy consumption 
per successfully communicated information bit is essential in effective design and deployment 
of wireless systems. This accounts for example for cases such as low-energy sensor networks 
in ecological environment monitoring as well as energy consumption in infrastructure devices 
of cellular systems. In addition, it is in line with global policies and strategies on low energy 
consumption and awareness on environmental issues which, among others, has led to the rapid 
emerge of green communications 1131, Pfl. 

It has been shown that multi-antenna systems constitute an effective method that can enhance 
spectral efficiency. However, this typically comes at a cost of complex transceiver circuitry and in 
massive systems with high energy consumption requirements. Furthermore, it is not currently fea¬ 
sible to embody large multi-antenna systems at hand held terminals due to spatial restrictions. As 
a result, cooperative communications have been proposed as an alternative solution that improves 
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coverage as well as performance under fading effects and have attracted significant attention due 
to their ability to overcome the limitations of resource constrained wireless access networks, 
see e.g. If5l- lfl8l and the references therein. A distinct feature of cooperative communications is 
that wireless agents share resources, instead of competing for them, which ultimately enhances 
system performance. In this context, various resource allocation algorithms and techniques have 
been proposed for improving the energy efficiency of resource constrained wireless networks. 
Specifically, the authors in lfl9l analyzed energy-efficient direct transmission adopting higher- 
level modulation for short distances, where circuit power is more dominant than transmission 
power. It was also suggested that high energy reduction can be achieved by optimizing the 
transmission time and the modulation parameters, particularly for short transmission distances. 
The authors in ff20ll addressed the optimal power allocation and throughput transmission strategy 
for minimizing the total energy consumption required to transmit a given number of bits. In 
mu, minimization of two-hop transmission energy with joint relay selection and power control 
was proposed for two policies: i) for minimizing the energy consumption per data packet; ii) 
for maximizing the network lifetime. In the same context, ||22 | - lf26j| addressed the modulation 
optimization for minimizing the total energy consumption for M —ary quadrature amplitude 
modulation (M— QAM), whereas energy efficient cooperative communication in clustered sensor 
networks was investigated in E71 . Energy efficiency in cooperative networks was also analyzed 
in E8l - lf35l by optimizing energy consumption based on the involved relay decoding strategy, 
modulation parameters, number of relay nodes and their distance from the source and the 
destination nodes. Likewise, an energy-efficient scheme was proposed in ll36l by exploiting the 
wireless broadcast nature and the node overhearing capability while an optimal energy efficient 
strategy based on the cooperative network parameters and transmission rate was reported in ll37l . 
Finally, the authors in |[38l analyzed realistic scenarios of energy efficient infrastructure-to-vehicle 
communications. 

It is also widely known that fading phenomena constitute a crucial factor of performance 
degradation in conventional and emerging wireless communication systems. Based on this, 
numerous investigations have addressed the effect of different types of fading conditions on the 
performance of cooperative communications If39l - ll56l . However, the vast majority of the reported 
investigations assume that the involved communication paths are statistically independent to each- 
other. Nevertheless, this assumption is rather simplistic as in realistic cooperative communication 
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scenarios the wireless channels may be spatially correlated, which should be taken into account 
particularly for deployments relating to low-energy consumption requirements. Based on this, the 
authors in 11571 addressed the spatial correlation in relay communications over fading channels 
whereas the performance of a decode-and-forward (DF) system with M— PSK modulated signals 
in triple correlated branches over Nakagami —m fading channels using selection combining was 
investigated in |[58l . In the same context, the authors in lf59l analyzed the performance of a 
multiple-input-multiple-output (MIMO) DF system with orthogonal space time transmission over 
spatially correlated Nakagami —m fading channels for integer values of m. The performance of 
a two hop amplify and forward (AF) relay network with beamforming and spatial correlation 
for the case that the source and destination are equipped with multiple antennas while the 
relay is equipped with a single antenna was investigated in |[60l . Likewise, spatial correlation in 
the context of indoor office environments and multi-antenna AF relaying with keyhole effects 
was analyzed in iffTHl and lf62l . respectively whereas the effects of spatial correlation on the 
performance of similar relaying systems were analyzed in |[63l - lf66ll . However, to the best of the 
authors knowledge, a comprehensive exact and asymptotic error rate analysis for regenerative 
systems over spatially correlated channels using maximum-ratio-combing (MRC) as well as a 
detailed energy-efficiency analysis and optimization, have not been reported in the open technical 
literature. 

Motivated by the above, the aim of this work is twofold: we, firstly, derive exact analytic 
expressions for the SER of a two-hop DF relay system over spatially correlated Nakagami —m 
fading channels for both M —QAM and M —PSK constellations along with simple and accurate 
asymptotic expressions for high signal-to-noise ratio (SNR) values. Secondly, we provide a 
comprehensive analysis of energy-efficiency and the corresponding optimization in terms of 
power allocation between cooperating devices. This is realized by minimizing the average total 
energy consumption of the DF relay network over multipath fading conditions, for a given 
destination bit error rate (BER) and maximum transmit power constraints. 

In more details, the technical contributions of the present article are outlined below: 

• Exact closed-form expressions are derived for the end-to-end SER of M —QAM and M —PSK 
based dual-hop regenerative relay networks with MRC reception at the destination over 
Nakagami —m multipath fading channels with arbitrary spatial correlation between source- 
destination (S-D) and relay-destination (R-D) links. 
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• Simple closed-form asymptotic expressions are derived for the above scenarios for high 
SNR values. 

• The offered analytic results are employed in a comprehensive energy optimization analysis 
based on minimizing the average total energy consumption of the overall regenerative relay 
network under a given QoS target, in terms of bit-error-rate (BER), and maximum transmit 
power constraints. 

The remainder of this paper is organized as follows: Section II presents the considered relay 
system and channel model while Sections III and IV are devoted to the derivation of the 
corresponding exact and asymptotic error rate results. The total power consumption models 
are presented in Section V while the analysis of energy minimization and power allocation 
optimization based on the given constraints are provided in Section VI. Section VII presents 
the corresponding numerical results along with extensive analysis and discussions while closing 
remarks are provided in Section VIII. 


II. System and Channel Model 


We consider a two-hop cooperative radio access system model consisting of source node (S), a 
relay node (R) and a destination node (D), where each node is equipped with a single antenna, as 
illustrated in Fig. 1. Without loss of generality, the system can represent both a conventional and 
emerging communication scenarios such as, for example, a mobile ad-hoc network or a vehicle- 
to-vehicle communication system. The cooperative strategy is based on a half-duplex DF relaying 
where transmission is performed using time division multiplexing. It is also assumed that the 
destination is equipped with MRC reception and that information signals are subject to multipath 
fading conditions that follow the Nakagami —m distribution!!. 

In phase I, the source broadcasts the signal to both destination and relay nodes and the 
corresponding received signals can be expressed as 




vs,d = 


p, 


Ls,D 


-h s ,DX + n s ,D 


( 1 ) 


and 


1 It is noted that the considered system requires the least resources in terms of bandwidth and power compared to multi-relay 
assisted transmission and thus, it can be adequate for low complexity and low-energy wireless networks. 
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Fig. 1. A dual-hop cooperative single relay model. 


ys,n — 



hs,RX + n s ,R 


( 2 ) 


respectively, where Ps is the transmit power, x is the transmitted symbol with normalized unit 
energy in the first transmission phase and Pls* and Pl b * denote the path loss values in the S-D 
and source-relay (S-R) paths, respectively. Also, h R ,n and h s R are the complex fading coefficients 
of the S-D and S-R wireless links, respectively, whereas risjj and ris.n are the corresponding 
complex Gaussian noise terms with zero mean and variance N 0 . The relay then checks whether 
the received signal can be decoded correctly, which can be, for example, realized by examining 
the included cyclic-redundancy-check (CRC) digits or the received SNR levels |[55l . If571 . Based 
on this, if the signal is successfully decoded, the relay forwards it to the destination during phase 
II with power P R = P R ; otherwise, the relay does not transmit and remains idle with P R = 0. 
Hence, the signal at the destination during phase II can be represented as 


Vr,d 



h Rt D% + n R: D 


(3) 


where P R is the transmit power of the relay, Pl rd is the path loss of the R-D path and h R D 
and n R j) denote the channel coefficient and complex Gaussian noise term with zero mean and 
variance N 0 , respectively. Finally, the destination combines the received direct and relayed signals 
based on MRC principle where the combined SNR can be expressed as follows l(57l . If67l 


7mrc — 


( Ps/Pl s , d ) I h s ,p +(Pr/ Pl R jD ) I h R)D 
No 


(4) 


The fading between the devices is assumed to follow the Nakagami —m distribution, which is a 
widely used model as the fading parameter m can easily account for both severe and moderate 








7 


12 i 12 12 

fading conditions. Thus, the corresponding channel power gains | h S) D \ ,| hs,R \ and | h RyD \ 
follow the gamma distribution lf 68 l with different power parameters, n 777 ’Ds~R’ 07777 - an d fading 
parameters m S} D, m s.n and rn RJ) , respectively. 

In the considered regenerative system, arbitrary spatial correlation is assumed to exist between 
the S-D and R-D paths, as also adopted in the semi-analytical contribution of lf57l . To this effect, 
the corresponding MGF for the case of Nakagami —m fading is expressed as |[69l 


M(s) = 1 - 


(1s,dI p Ls, d f 7 r,d/Pl r ,d) . (1 P)1s,d1r,d/{ p l s , d Pl R:D ) 


-s + 


m 


m~ 


s < 0 

(5) 


where D and 7 R D are the corresponding average SNR values, s denotes the MGF parameter 
and p = Cov( h s .n , \h R:D | 2 )/V Var (l h s, D | 2 )Var(|/i 7?iD | 2 ) represents the involved correlation 
coefficient ll 68 l . If69l . with Cov(-) and Var(-) denoting covariance and variance operations, 
respectively. It is acknowledged that other practical impairments, such as co-channel interference, 
are not considered in this work. 


III. SER for M- QAM Modulation in Nakagami-tti Fading with Spatial 

Correlation 

This section is devoted to the error probability analysis of the dual-hop cooperative network 
over Nakagami —m fading channels with spatial correlation between the direct and the relay- 
destination paths ll 68 l . To this effect and assuming MRC reception, the average end-to-end SER 
can be expressed as If57l 


zC 


SERn =F ( 


D r QAM 


( P S / P L S d )^S,D 9QAm\ 1 


1 + 


Nom c sin 2 9 J 


QAM 


(P’s/ PL SR )^S,R9QAM\ mS ' R 


1 + 


Aq ms tR sin 2 9 J 


F, 


QAM 


( \ 1 
[(Ps^S.d/Pl s d )+(Pr^R,d/ Pl r d )\ 9QAM (1- p) P SP r^S,D^R,D9q AM \ 

N 0 m c sin 2 9 N o p L s r> Pl r £> m c sin4 9 J 


x <1 — F ( 


QAM 


/7 (P’s/ P L S}R ) tt S ,R 9QAM \ mS ’ R 

\ Noms lR sin 2 9 J 


(6) 


























where m c = m s , D = m R D , g qAM = 3/2 (M - 1) and 

4 / i \ f 71 / 2 4 / l \ 2 W 4 

^ Ml vm = -(1-^)1 l m*. m 

The first two terms in © refer to the cases of incorrect and correct decoding of the received 
signal at the relay node, respectively, whereas the integral representation in © is used for 

_ (J 

evaluating the SER D numerically. In what follows, we firstly derive a closed-form expression 
for the average SER in the case of direct communication mode. This expression is subsequently 
employed in the derivation of exact closed-form expressions for the average SER of M—QAM 
and M— PSK modulated regenerative systems over Nakagami— m fading channels with spatial 
correlation. Furthermore, it is used in the analysis of the energy consumption model and energy 
minimization in Sec. VI as it allows the derivation of an accurate expression for the energy 
consumption in the direct transmission, which acts as a benchmark in the evaluation of the 
energy reduction of the cooperative system. 


A. Exact SER for the Direct Transmission (DT) 

Theorem 1. For P R , Pl s d , &s,d, No, #qam £ M + , Me N and ms,D > the symbol error rate 
of a M—QAM direct transmission scheme can be expressed as follows, 


-D 


ser d = 


y/nM(ms,D N 0 P Ls d + Ps^s,DgQAM) ms ' D 
r (ms d + 4) ( 1 

x < TV ’ I 2 Fl m S,Di m S ,D + 1, 

1 (m +1) V 2 

\[2(\[M — 1) 


m s ,DN 0 P Lg 


mS,DN 0 PL StD + Ps&S,DgQAM 


71 




3 1 


m s ,DN 0 P Lsj 


2 2 2(ms,DNoPL s , D + -Ps , ^5,d5 , qam) , 


where 2 -Ei(.) and F\ (.) denote the Gauss hypergeometric function and the Appell hypergeometric 
function of the first kind, respectively. 


Proof: The proof is provided in Appendix A. 











B. Exact SER for the Cooperative-Transmission (CT) 


In this subsection, we derive a novel closed-form expression for the average SER of the co¬ 
operative transmission scenario when the involved relay node decodes and forwards successfully 
decoded information signals to the destination. To this end, it is essential to firstly derive exact 
closed-form expressions for two important indefinite trigonometric integrals. 

Lemma 1. For a, b, rri E M + and 2m — | E N, the following closed-form expression is valid, 


J(a, b,m) = 


rd£ 


1 -I- a 4- b ) 

1 ^ sin 2 (0) ^ sin 4 (0) J 

(2m - ^ ^ (a + 2 sin 2 (6 l ) - -fa 2 - 4 b) (a + 2sin 2 (0) + -J a 2 — 46) ’ 
V ^ / (2 + a — -fa 2 — 46) m (2 + a + -fa 2 — 46) m 

, T7 (7 1 1 7 | 3 2 cos 2 ( 0 ) 2 cos 2 ( 0 ) A 

COS 1+2l (9) y + 2 ,m,m,l + ^, 2+a _p a S_ 4b , 2 +a+Va 2 -4b) 


X 


(1 + 2/) (sm 4 (8) + asm 2 (8) + b) 

Proof: The proof is provided in Appendix B. 


+ C. 


(9) 


Lemma 2. For a, b, rri, n G M + and rri + n — 1 G N, the following closed-form expression is 
valid 


/C(a, b, m, n ) = 


1 + 


sin 2 (0) 


1 + 


b 

sin 2 (0) 


r d 8 


1=0 


m + n 
l 


_ in (-1)' cos 1+2i (») F, (i + i,m, n, ( + |, ) 


(1 + 2/)(l + a) m (l + b) r 


+ c. 


( 10 ) 


Proof: The proof is provided in Appendix C. ■ 

To the best of the authors’ knowledge, the generic solutions in the above Lemmas have not 
been previously reported in the open technical literature. These results are employed in the 
subsequent analysis. 


Theorem 2. For {Pg, Pr, Pl S}D , Pl s , r ,Pl R}D Ms,d,^s,rMr,d, N 0 } G M + , MgM, tos,d > \, 
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mg,R > \, m R,D > ms,D, 2 m c — ~ G N and 0 < p < 1, t/ie S£R of M—QAM based DF 
relaying over spatially correlated Nakagami—m fading channels, can be expressed as follows: 


rC 


SER^ = 


2 (m c -i)! 2 Fi ^ - m c , m c , 2 +kr) 


x 


X 


7im c \ M{pfM — 1 ) 1 (1 + ai) mc ttM(VM-I) 2 (1 + ai) mc 

2 (pT>S,R — 2) • 2-^1 (jtlS,Rj 2’ m S,R T 1 ? 1^7) 2 y r 2 F\ ^2 > 2 m S,Ri m S,Ri 2’ 2’ 2+2bi ) 


+ 


+ < 1 


y/nm StR \ MQ/M - I)" 1 (1 + 7 r M (\/M - l) -2 (1 + 6i) ms '« 

2 - |)! 2 Ei (^rris.Ri \,m s , R + 1 , 2 ^^! (§, § - m s , R , m s , R , §, §, 27 ^ 7 ) 


7T ,r! M{y/~M — 1) 1 (1 + bf) ms < R 


7T M (VM - l)- 2 (1 + 6i) m s.fl 


^ / 2m e - (~1) ; 4(7 m - l) 2 2 m -- ; -^F 1 (l + 1, m c , m c , / + f, 2.A, 2 B) 


1 =0 

2m c — 


Y' f 2 m c - |\ (-1)'4(a/M - l)Ei (/ + l,m c , m c , 1 + %,A, B) 
^ \ l ) tiM (1 + 2l)(a 1 d 1 ) m -(l - A)~ mc (l - B)~ m c 


«=o 


( 11 ) 


where a\ — Ps^s,d9qam / {PLs, D m s,DNo), hi — Ps^s,r9qam / (PL S:R Noms,R)> c i — Pr^r,d9qam/ 
( PL R}D No m R,D ), G?1 = (1 ~ p)PR^R,DgQAM/{PL R ^ D No m R,D) and 


0 = 


( 12 ) 


2 + cii + C\ {-)_}• \J (01 + C 1) 2 — Aaidi 

Proof: The first term in © corresponds to the direct transmission and thus, it can be 
expressed in closed-form based on Theorem 1. Likewise, the second and the fourth term in © 
have the same algebraic representation as ([62b and (l63l) in Appendix A. Therefore, they can be 
readily expressed in closed-form by making the necessary change of variables and substituting in 
(1681) and (T7Ql) . As for the third term in ©, it is noticed that it has the same algebraic representation 
with ©. As a result, a closed-form expression is deduced by determining the following specific 
cases in ©, which practically evaluate ©, 


J (a,h,m, 0, {^ 4 }) 


■{#} 


1 + 


sin 2 (e») 


+ 


b 

sin 4 (0) 


rd e. 


(13) 






















Therefore, by carrying out some long but basic algebraic manipulations and substituting in ([6]) 
along with the aforementioned closed-form expressions, one obtains (flTT) . which completes the 
proof. ■ 


Remark 1. Equation ( fiTI ) reduces to the uncorrelated scenario by setting p = 0. However, an 
alternative expression for this case which is valid for the case that {p + m R .p — 1 } G N 
can also be deduced by applying the derived expressions in Theorem 1 and Lemma 2 in / l5Zl 
eq. (11)], namely, 


c 


SER^ = 


2T (m s ,D + |) 2+1 (m s ,D, §, 1 + m s , D , 2y/2(F 1 (j, | - m s , D , |, 


x 


+ < 1 


V^M(VM-l) ^(1 + m S) D){f + ai) ms ’ 

f 2 r (m s ,R + 2) 2+1 2,1 + m s ,R, ,+ fel 

D 7 t(v/M - 1)- 2 M(1 + ai) ms . D 

) 2 \J~ 2 F\ ^ 2 , 2 m S,Ri 2 ,2’ 2+261) | 

v /?M(v / M - 1) _1 T(1 + ms, R ){l + 

2 T (m s ,R + 2) 2+1 1 + m s ,R, 

: tt(\/M - 1)- 2 M(1 + 

2v / 2+i (|, 5 m StR , 2,|, 2+261) 
+--- 



X 


y/HM(sfM — 1) 1 r(l + m s ,/ ? )(l + 6i) ms ’« 7 t(\/M — 1) 2 M(1 + bi) ms > R 

i 4(-l)'(+M - 1)F, (i + i , ms, D , '"Hu,, l + tE) 


E 

1=0 


+ ai) m ^(l + ci) m ^(l + 21) ( m s , D + m R , D + \) 


+ 


m S , D +rn«. D -i (_1 )Z 2 H( VM - l) 2 Fi (/ + m 5 ) D, m R,Di l + |, 


1 


1 


E 

Z=0 


2+2ai ’ 2+2ci 


l\nM(l + ai) m ^(l + ci) m «^(l + 21) (m s , D + rftR.D + ±) 


(14) 


where ( x) n = T(x + n)/T(x) denotes the Pochhammer symbol fi73\l . 


C. Asymptotic SER for the Cooperative-Transmission 

Simple asymptotic expressions can be derived for the case of high SNR at the three paths 
of the system. To this end, it is essential to firstly derive a closed-form expression for another 
trigonometric integral. 

Lemma 3. For rri g E, the following generic closed-form expression holds, 
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Proof: The proof is provided in Appendix D. ■ 

Lemma 3 is subsequently employed in the derivation of the following proposition. 

Proposition 1. For {Ps, Pr, Pl s ,di Pl S ' R , Pl R ' D ,Qs,d,&,s,r,&,r,d, N 0 } e M + , Me N, ms,D > 
> 4, > 4, 2m c — 4 G N and 0 < p < 1, f/re SLR of M—QAM based DF 

relaying over spatially correlated Nakagami—m fading channels in the high SNR regime can be 
expressed as 


_ ( P Ls ,oN 0 m c \ mc (, _ J_\ f 2r ( m c + |) J_\ 2^1 (4,l,m c + f,-l) | 

D \Ps^s,d9qam) V VMj | v^r(l + m c ) \ VM) vr2 m =" 2 (1 + 2m c ) J 

x / _ J_\ f 2r (m g|fl + j) _ / _ J_\ 2 *\ & + j,-l) | 

V Ps^s,r9qam) v vW \V^r(l + m 5)R ) \ v/M/ n2 m s,R~ 2 (l + 2m s>R ) J 


P^S.d Plr,d ^0 "»c 


(1 - p ) P s PrFIs,d^r,d9qam 


X 


1 - 


JJ) 

v/f; 


\ 2 ^ ( 2m c + 2 ) | 

fi 1 ^ 

2 F 1 (i,l, 2 m c +|,-l)'l 

[ vFr(i + 2 m c ) 


7T 2 2rric - 2 (l + 4m c ) J 


( 16 ) 


Proof: In the high SNR regime, it is realistic to assume that Ps^s. D >> PL StD No, Ps^S,R » 
P Lsr N 0 and P R VL R)D » P Lrd N 0 . Based on this, an integral representation was formulated in 
m eq. (28)], 


r C 


SER d ~ A r A 


S,R 


N 0 m c 


■Y 


Ps^s,d9qam J \Ps^s,r9qam J 


N 0 m s . 


R 


-) 


ms,R 


+21 


2c 


iv 0 v 


(1 — p)PsP r FIs,dQr,d9qam 

(17)' 


where 



4 / i \ W 2 4 / I \ 2 r + 4 

A gfl = - 1 - —= / sin 2ms ’« d0-( 1 — -f= / sin 2ms ’« d0. (19) 

’ ^ V Vm) Jo vt V VmJ Jo 
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Evidently, the terms A c , A 2c and A S p can be expressed in closed-form with the aid of Lemma 
3. Based on this, by performing the necessary change of variables in and substituting in 
(fTTT) . equation (fT6l) is deduced, which completes the proof. ■ 

Remark 2. Using ( 1 1 61) . the correlation coefficient for the case of M— QAM modulation can be 
expressed in terms of the corresponding source and relay powers, fading parameters and average 
SER, namely 


P 


/ 

$ER° d 

V 


I<:>, 


2CT(2m c +i) 
V / 7rr(l+2 m c ) 


C 2 2 Fi(i,l,2m c +| ,-l) 
7r 2 2m c _2 (l+4m c ) 


l 

m c 


j 2CT ^ m c + 

C 2 2^1 ( ^ >l> m c + § , — l) 

1 j 2Cr^m SjJ j + ^) + 1 

| \/^ r (l+ m c) 

7,2 m c -‘2 (l-(-2 me) 

J [ + n 2 rn S , R - 2 ( l + 2 m S R ) j 


K^ c K 2 




1 

m c 


( 20 ) 


where g = g QA M,C = (l - 1/y/Htj ,K\ = (N 0 m c P Ls D /P s CL s , D g ), K 2 = (N 0 m S:R P LStR /P s ^s,r 9), 
and K 3 = (N^m 2 c P LsD P LRD /PsP R n s ,D^ R ,Dg 2 )- 


IV. SER for M-PSK Modulation in Nakagami-tti Fading with Spatial 

Correlation 

Having derived novel analytic expressions for the case of M— QAM modulation, this section 
is devoted to the derivation of exact and asymptotic closed-form expressions for the case of 
M— PSK constellations. 

A. Exact SER for the Cooperative-Transmission 

Theorem 3. For {P R , Pr, Pl S iD , Pl s , R : Pl RiD ,^s,d^s,r,^r,d, No} 6 M + , M e N, rnsjj > 
ms,R f §> m n,D > \, 2 m c - | 6 N and 0 < p < 1, the SER of M—PSK based DF 

relaying over spatially correlated Nakagami—m fading channels, can be expressed as follows: 
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where <22 — Ps^s,DgpsK/{PL StD fns,DNf), b 2 — Ps^s,r9psk/(PL s , R Noms, R ), c 2 — P r Q R) d9psk / 
(PL RtD No'm R ,D), d 2 = (1 — p)P R Q R ,DgpsK/{PL R ' D N 0 m RtD )- 


Proof: As a starting point, the average SER of M— PSK modulated DF systems over 
Nakagami—m fading channels with spatial correlation can be formulated as follows lf57l eq. 
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The above four integrals have the same algebraic form as the integrals in Theorem 1 and Lemma 
1. Thus, the proof follows by performing the same necessary change of variables and substituting 


in 


B. Asymptotic SER for the Cooperative-Transmission 

Proposition 2. For {P s , P R , P LStD , Pl S:R , Pl RiD Ms,dMs,r,^r,d, No} eR + , M e N, m s , D > 
2, > §, m R,D > \, ms,D, 2 m c — | G N and 0 < p < 1, the SER of M—PSK based DF 

relaying over spatially correlated Nakagami—m fading channels in the high SNR regime can be 
expressed as 
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Proof: The asymptotic SER for high SNR values was formulated in Il57l eq. (27)] 
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Notably, the integrals in (1261 ) and (1271) have the same algebraic representation as the integral in 
Lemma 3. As a result, by performing the necessary change of variables and substituting in (1251) . 
one obtains (l24l) . which completes the proof. ■ 

Remark 3. Based on (l24l) . f/ze corresponding correlation coefficient can be expressed in terms 
of the corresponding source and relay powers, fading parameters and average SER as 


P 


= 1 - 


( ( r(m e + ^) col^) 2 F l(^-?- m e-|. cos2 (in-)) | f r ( m S,fl+^) cos (l w) 2 F 1 ( 2 ' 3 ' §- cob2 ( TT ) ) ) \ mc 

| 2-^/tF m c ! ' 7r [ | 2\pR mg tt 


SER 




(28) 




where g = gpsK is set in the I\\, I\ 2 and K : > terms, which are given in Remark 2. 


Figure 2 illustrates the SER performance as a function of SNR for 4—QAM/QPSK modula¬ 
tions. The source-destination transmission distance is indicatively considered at 600m while the 
relay is assumed to be located in the middle and the transmit power is shared equally to the 
source and the relay. The corresponding path-loss effects are considered by adopting the path 
loss (PL) model in f70i1 namely 


PL i}j [dB\ = 148 + 40 log 10 (dij [km]) 


(29) 
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which has been shown to characterize adequately harsh communication scenarios and is partic¬ 
ularly applicable to mobile relaying and device-to-device communications. It is clearly observed 
that the empirical simulated results are in excellent agreement with the respective analytical 
results. Furthermore, the simple asymptotic results are also highly accurate at higher SNRs. 



Fig. 2. Example SER performance over Nakagami— m fading channels with ms,D = ms,R = RIr,d = m = 
{0.75, 1.25}, = tls ,r = Qr,d = OdB for 4—QAM/QPSK constellations and different values of spatial correlation. 

V. System Power Consumption Model and Analysis 

In this section, motivated by the general interests towards green communications and increas¬ 
ing incentives to save energy, we quantify the total energy consumption required to transmit 
information from the source to the destination. We assume that the transceiver circuitry operates 
on multi-mode basis i.e. i) when there is a signal to transmit, the circuits are in active mode; ii) 
when there is no signal to transmit, the circuits operate on a sleep mode; in) the circuits are in 
transient mode during the switching process from sleep mode to active mode. The elementary 
block diagrams of the assumed transmitter and receiver are illustrated in Fig. 2 and Fig. 3, 
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IQ Mixer RF PA \ / 



Fig. 3. Elementary direct-conversion transmitter. 

respectively. This model is based on the energy and layout area efficient direct-conversion 
architecture which is commonly used in wireless transceivers. It is also assumed that all nodes 
are equipped with similar transmitter and receiver circuit blocks and that the power consumption 
of the active filters at the transmitter and receiver is similar. 

Considering a node that transmits L bits and total transmission period T, the transient duration 
from active mode to sleep mode is short enough to be neglected. However, the start-up process 
from sleep mode to active mode may be slower due to the finite phase-locked loop (PLL) 
settling time in the frequency synthesizer. By denoting the duration of the sleep, transient and 
active modes as T sp , T tr and T on , respectively, the total transmission period is defined as T = 
T sp + T tr + T on , with T tr being equal to the frequency synthesizer settling time. Based on this, 
the total energy required to transmit and receive L information bits is expressed as 

E = P on Ton + PspTgp + PtrTtr (30) 

where P on , P sp and P tr denote the power consumption values during the active, sleep and transient 
modes, respectively. In realistic circuit designs, the power consumption in the sleep mode can be 
considered negligible compared to the active mode power |[T9ll and thus, P sp ~ 0. It is also noted 
that power consumption during the transient mode practically refers to the power consumption 
of the frequency synthesizers. Based on this, it is assumed that P tr = 2 1’to ; therefore, using 
the power consumption values at both transmitter and receiver sides during the active mode one 
obtains 
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Fig. 4. Elementary direct-conversion receiver. 


Pan = P< 


ont 


P 

1 n r i 


(31) 


where P ont is the total transmitter power consumption that accounts for the sum of signal 
transmission and transmitter circuit powers and P onr is the total receiver power consumption. 
Hence, it follows that, 


Pont ~ Pt + Pamp. + PcT x (32) 

and 

Ponr = PcR x (33) 

where P t is the signal transmission power, P a m P is the power consumption of the RF power 
amplifier and Pct x and Pcr x denote the total transmitter circuit power and the total receiver 
circuit powers, namely, 


Pcr x — PdSPtx + PdAC + PfU + PMix + PLO 


(34) 


and 


PcR x = PdSP Rx + PadC + PvGA + 2 Pfu + Pmix + P LO + Plna (35) 

respectively. The P C t x measure consists of the following power consumption entities: digital 
signal processor (DSP), Pdsp t P digital to analog converter (DAC), Pdac ; active filter, Pfu', 
IQ Mixer, Putx and synthesizer, Ppo- Likewise, the active power consumption at the receiver 
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comprises the power consumption values for digital signal processor (DSP), Pdsp R x \ analog 
to digital converted ADC), Padc', variable gain amplifier (VGA), Pvga ; active filter, P Fii ; IQ 
Mixer, Pmix\ synthesizer, P LO ', low noise amplifier (LNA), Plna If75l - Based on this, the total 
required circuit power consumption is given by 


Ptc — Pct x + Pcr x ■ (36) 

It is also noted that for signal transmission power P t , the power consumption of the RF-power 
amplifier can be modeled by P amp = aP t , where a — | — 1, with // and £ denoting the respective 
drain efficiency of the amplifier and the peak-to-average power ratio (PAPR), which depends 
on the modulation order and the associated constellation size. Based on this, for the case of 
square uncoded M —QAM modulation, £ = 3 ^j|~* and T on = Pp = -P, where b = log 2 M is 
the constellation size, L is the transmission block length in bits and T s is the symbol duration 
which relates to the bandwidth B as T s « T lf23l 


VI. Energy Optimization and Power Allocation 

In this section, we deploy and combine the results of the previous sections and analyze the 
total energy required to transmit information efficiently from the source to the destination. To 
this end, we firstly quantify the total energy consumption in the direct communication scenario. 
Hence, by applying ([30b . (l32l) and (03b and recalling that P sp ~ 0 and P tr = 2 P LO , the average 
energy consumption per information bit is given by |[76ll 

_ -pD _ ((1 + a)P S + PcT x + PcRx) Ton + 2 PloTit 

h/rp h/rp ^ W / ) 

where Ps denotes the source transmit power. In order to determine the average total energy 
consumption in the corresponding cooperative transmission system deploying the DF protocol, 
we formulate the total average power consumption, which is a discrete random variable that can 
be statistically expressed as 


pC _ 
r T — 


Pct x + (1 + a) P s + 2P C r x , 
Pct x + (1 + a) P R + P C r x , 


with Pr = 1 


(38) 


with Pr = 1 — SER 


■S,R 
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where Pr denotes the relay transmit power. The first term of (1381) refers to the absolute total 
power consumption by the nodes in the first transmission phase, while the second term represents 
the power consumption in the second phase, subject to correct decoding of the received signal 
by the relay, which is indicated by the probabilistic term (l — SER^r). Hence, the average total 
power consumption in the cooperative transmission mode can be expressed as 


Pt — Pcr x + (1 + Ot)Ps + 2Pcr x + ( Pct x + (1 + (x)Pr + Pcr x ) (l — SER^r) . (39) 

Based on this, the corresponding average energy consumption per information bit is given by 

E$ = PtTou + 2PL ° Tt R (40) 

The achieved energy efficiency enhancement by the cooperative transmission is determined 
with the aid of the cooperation gain (CG), which is the ratio of the energy efficiency of 
cooperative transmission over the energy efficiency of the direct transmission, per successfully 
delivered bit, namely, 


CG = 


e d t (1 


BER 


E% (1 


-D 


(41) 


BERr, 

Evidently, when the resulting ratio is smaller than one, it indicates that direct transmission is 
more energy efficient and thus, the extra energy consumption induced by cooperation outweighs 
its gain in decreasing the BER of the system. 

In what follows, the above expressions are employed in formulating and solving the energy 
optimization problems aiming to guarantee certain QoS requirements, namely, target destination 
BER. In the same context, we additionally provide the optimal power allocation formulation for 
the cooperative transmission scenario under the maximum total transmit power constraint. 


A. Direct Transmission 

We first consider the energy optimization problem for minimizing the average total energy 
consumption in the direct communication scenario with the maximum transmission power and 
target bit error rate, p*, as constraints. We assume that the power consumption of the circuit 
components are fixed and independent of the optimization. Thus, the only variable in the 
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optimization is the transmit power of the source. To this effect and with the aid of ([37l) . the 
optimization problem for the direct transmission mode can be formulated as follows: 


min E T 

Ps 


(42) 


subject to: P s < P maxt: P s > 0 


BER,£ =p*. 


Deriving the minimum average total energy required in the direct communication scenario, 
requires prior computation of the corresponding symbol error probability. This is realized with 


the aid of ® which is expressed in closed-form in terms of 2 Fi yn, f; m + 1; iq—J and 
Fi 2 — m, m, §; 2+ |^ j function!]. It is recalled that these functions are widely employed 

in natural sciences and engineering and their computational implementation is rather straight¬ 
forward as they are built-in functions in popular software packages such as MATLAB, MAPLE 
and MATHEMATICA. It is also noted that the representation of these functions in the present 
analysis allows the following useful approximative expressions: 2 F\ (m, m + 1; yjyy'j — 1 and 
F\ 2 — m, m, |; 2+ \ ai j ~ F\ 2 — m, m, |; 0). The accuracy of these approximations 

is validated through extensive numerical and simulation results which indicate their tightness for 
random values of m and moderate and large values of a\. To this effect, the following accurate 
closed-form BER approximation for M—QAM modulated signals is deduced 



Importantly, the Appell function in (l43l) can be expressed in terms of the Gauss hypergeometric 
function, 



(44) 


As a result, equation (l43l) becomes 



2 


(45) 


2 For the sake of simplicity, we assume that m = ms,D■ 
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It is evident that (1451) is a function of the modulation order, the severity of multipath fading and a i . 
Therefore, by substituting the targeted QoS p* in (l45l) . recalling that Gq = (Ps^s,d9qam) /(A t qPl sd ) 
and carrying out some algebraic manipulations, one obtains 


where 


Po ~ 


™>n q p LS iD 

^s,d9qam 



(46) 


C = 


4(/M- 1) 
M log 2 M 


(VM T(m + |) 


(47) 


a /2 7T 2^mT(m) 

To this effect and with the aid of (1371) and (l46l) . it follows that the minimum total energy per 
information bit required for direct transmission for meeting the required QoS can be expressed 
in closed-form as 


—d* _ ( Pct x + Pcr x ) T on : (1 + a)N 0 ms t DT on PL SiD 

f J rj-\ - | 

L L Qs,d9qam 

Based on the total energy consumption in (1481) and given the constellation size b = it 

is shown that the proposed energy expression comprises the transmission energy E t and circuit 
energy E c , namely, 


1 /m 
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+ EELl, ( 48 ) 
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where C can be expressed as a function of the transmission time T on as follows 
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Hence, by inserting (l50l) in (|49|) . the following analytic expression for the transmission energy 


per information bit is deduced 
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Likewise, the circuit energy, E c , can be expressed as 

Ec = ( PcT x + PcR x ) Ton ■ (52) 

Notably, equations (l50l) and (1511) indicate that for a fixed bandwidth B and packet length L, the 
transmission energy is a decreasing function with respect to the product T on B whereas the circuit 
energy increases monotonically with respect to T on . In addition, it is shown that the transmission 
energy is dependent upon the transmission distance ds,D and the severity of fading m, while the 
corresponding circuit energy remains fixed regardless of the value of ds. n and m. 

B. Cooperative Transmission 

This subsection presents the energy optimization and power allocation problem when the 
involved relay forwards successfully decoded signals, generally at different power than the power 
of the source. Evidently, the respective optimization model is a two dimensional problem and 
thus, we formulate the energy minimization problem with two optimization variables, namely, 
the source transmit power, P$, and the relay transmit power, P R . In this context, the aim is to 
minimize the total energy consumption of the overall network instead of minimizing the energy 
consumption at individual nodes. Based on this and with the aid of (l40l) . the optimization problem 
can be formulated as follows: 


min E R (Pg, P r ) 

Ps,Pr 

subject to: (P s + Pr) < Pmaxt, Ps > 0, Pr > 0 ( 53 ) 

mR C D (P s ,P R )=p*. 

The above optimization task is a non-linear programming problem since the objective function 
and the constraint BER are both non-linear functions of P s and P R . It is also recalled that 
Karush Kuhn Tucker (KKT) conditions that handle both equality and inequality constraints are 
in general the first order sufficient and necessary conditions for optimum solutions in non-linear 
optimization problems (NLP) provided that certain regularity conditions are satisfied. To this end, 
using the Lagrange multipliers Ai and A 2 , A 3 and A 4 , for the equality and inequality constraints 
we set the corresponding Lagrangian equation that depends on the optimization variables and 
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multipliers while meeting the KKT conditions in IfTTI for the non-linear convex optimization 
problem, namely 

L = Et + \i (bER^ - p*) - A 2 P S - A 3 P R + A 4 {{Ps + Pr) - Pmaxt). (54) 

The proof for the convexity of the optimization problem is provided in Appendix E. 

Based on (l54l) . the KKT conditions for the problem can be expressed as follows: 

WErp + AiVBER d — A 2 V-Ps 1 — A 3 V Pn + A 4 V (Ps + Pr) — 0 (55) 

whereas the associated complementary conditions are given by 


BERp — p*, Ps + Pr < Pmaxt 

Ai (bER^ - p*^j = 0 , A 2 Ps = 0, A 3 P R = 0 and A 4 (P s + Pr- P max t ) = 0 (56) 

Ai, A 2 , A 3 , A 4 > 0. 

In the above set of complementary KKT conditions both A 2 and A 3 represent inactive constraints 
and therefore, they can be assumed zero. To this effect, by applying (1551) and setting the 
derivatives w.r.t Ps and Pr to zero, the following useful set of equations is deduced 


and 
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dP s 
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r c 


dE T 9BER d 

+ Ai D + A 4 = 0. 


(58) 


OPr ()Pr 

Solving for A 4 from (l57l) and substituting in (l58l) yields the following relationship which depends 
only on one of the Lagrangian multipliers 
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Based on this and using the fact that Ai > 0, one obtains the following necessary condition for 
minimizing the total average energy consumption of the cooperative transmission mode at the 
optimal power values 

9E$(PS,Pi) ^ dE$(P s -,P£) 

m - w R ' (60) 

_ Q 

For a feasible set of optimal powers, the BER D = p* and P s + P R < P maxt constraints must be 
satisfied. 

Analytic solution for the optimal powers in (1601) is intractable to derive in closed form. 
However, this can be alternatively realized with the aid of numerical optimization techniques, 
which can determine the optimal powers at the source and relay nodes that minimize the 
average total energy consumption. To this end, we employ the MATLAB optimization tool 
box and its function fimincon in the respective numerical calculations for allocating the available 
power optimally under the given constraints. Thus, the derived expressions and offered results 
provide tools to understand, quantify and analyze how much energy can in general be saved, 
per successfully communicated bit in the system, if transmit power allocation and optimization 
beyond classical equal power allocation is pursued in the cooperative system, on one side, and 
how much energy can be saved against the classical non-cooperative (direct transmission) system, 
on the other side. Furthermore, the considered values in the present paper are indicative and are 
selected in the context of demonstrating the validity of the proposed method. Therefore, the 
derived optimization flow can be readily extended to arbitrary design constraints for the total 
network power consumption and target destination error rate in the presence of Nakagami —m 
multipath fading conditions. 

VII. Numerical Results and Analysis 

This section demonstrates and evaluates the average total energy consumption of the considered 
regenerative system assuming that the S-D and S-R links are statistically independent whereas 
the S-D and R-D paths are spatially correlated. As a realistic example, we assume M —QAM 
modulation scheme over the S-D, S-R and R-D links, in case of cooperative transmission mode, 
and over the S-D link in the case of only direct communication. For the sake of simplicity, it is 
also assumed that all wireless channels are subject to Nakagami— m multipath fading conditions 
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x 10“ 6 



Fig. 5. Energy consumption per information bit versus source-destination distance for different relay locations over uncorrelated 
Nakagami—1.25 at target BER of 10 -2 for 4—QAM / QPSK constellation. 


with fls,D = -hsji = -hf.D = OdB. The involved path-loss effects are modeled by an example 
model of P Lij [dB] = 148 + 40 log[km]), which is also used in device-to-device based 
communications lITOll . and thus applies to mobile relaying as well. Furthermore, in order to 
simplify the geometry-related calculations, we assume that all nodes are located along a straight 
line, which satisfies the distance relationship d s . n = ds.R + d RD . However, it is recalled here 
that the path-loss and distance assumptions are only indicative in the context of the considered 
examples, while the provided analysis and optimization frameworks are valid more generally. In 
this context, we further assume the following system parameters: N 0 = — 174dBm/Hz; 77 = 0.35; 
Ttr — 5ns, L — 2 kbits, and P LO = 50mW ED, m. We also use the constant circuit powers 
as Pct x = lOOmW, Pcr x = 150mW and the maximum transmission power P max t = lOOOmW. 
The bandwidth of the system is assumed to be B = 200kHz and the noise figure Nf = OdB. 
Due to the linearity requirement of the M —QAM signals, the value of the drain efficiency is 
assumed 77 = 0.35, which is a practical value for class-A and AB RF power amplifiers (PA). 
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The considered system parameters are depicted in Table I and are used unless otherwise stated. 

TABLE i 

Assumed System Parameters 


N 0 = -174dBm/Hz 

Nf = 6dB 

B = 200kHz 

L = 2kbits 

Pct x — lOOmW 

Pcr x = 150mW 

P LO = 50mW 

Pmaxt = lOOOmW 

77 = 0.35 

Ttr = 5 /IS. 


We commence by analyzing the minimum energy per information bit required for the direct 
and cooperative transmissions when the relay node is taken into account and placed in different 
locations. The location of the relay node is represented with parameter / = ds,R/ds,D • Fig- 5 
illustrates the total energy consumption per information bit as a function of the transmission 
distance from source to destination for 4—QAM/QPSK with fading parameter of m — 1.25, 
destination target BER of 10~ 2 and zero spatial correlation under the maximum transmit power 
constraint. The transmit power allocation is carried out by the derived OPA scheme resulting 
to the indicative values in Table II. It is observed that distance thresholds separate the regions 
where DT performs better than CT and vice-versa. Furthermore, it is shown that when the 
relay is located in the middle, i.e. the source-relay distance equals the relay-destination distance 
(/ = 0.5), renders the best energy efficiency among all relay locations. This indicates that the 
configuration is almost symmetric in the source-relay and relay-destination distances, which 
assists the system to operate robustly in transmission over severe fading conditions. However, it 
is shown that at relatively small distances (here 0 < ds. d < 170m), the exact location of the relay 
does not affect substantially the performance of the cooperative system as it appears to remain 
almost the same in all considered scenarios. This renders the relay positioning and planning 
rather simple when the relay falls within this range, while it additionally provides insight e.g. 
for relay selection algorithms in the case of randomly distributed relays in emerging relay-based 
wireless networks. 

The corresponding energy efficiency is also analyzed for target BERs of 10~ 2 and 10 -3 . Fig. 6 
illustrates the average total energy consumption per information bit for 4—QAM/QPSK in both 
direct and cooperative transmissions in Nakagami—1.25 fading conditions under the maximum 
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Fig. 6. Energy consumption per information bit versus source-destination distance when the relay is located in the middle over 
spatially correlated Nakagami—1.25 for 4—QAM/QPSK constellations with different target BERs. 


transmit power constraint and the following spatial correlation scenarios: p = {0, 0.5, 0.9}. Also, 
the relay node is located in the middle and the transmit power is allocated optimally to the source 
and relay nodes in all cases. It is observed that for the fixed target BERs of 10' 2 and Hr 3 , 
the direct scheme outperforms the cooperative transmission only at average source-destination 
distances below 240m and 150m, respectively. On the contrary, for average distances greater 
than 240m and 150m, CT becomes more energy efficient as the transmit power constitutes a 
significant share of the average total energy consumption even under the worst spatial correlation 
scenario. Furthermore, it is shown that for the given target BERs the DT schemes attain maximum 
transmission distances of 390m and 250m, respectively, under the given maximum transmission 
power constraint while in both cases the cooperative transmission schemes extend to substantially 
longer distances. However, these advantages vary according to the level of the involved spatial 
correlation where the improvement in energy efficiency is inversely proportional to p, in both 
scenarios. The reason is that for every step of transmission distance, greater proportion of power 
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TABLE II 

Optimal Transmit Power Values for Source and Relay for Different Relay Locations over 
Uncorrelated Nakagami— 1.25 Links at target BER of 10~ 2 for 4—QAM / QPSK Modulation. 



Relay close to S (/ = 0.1) 

Relay in middle (/ = 0.5) 

Relay close to D (f = 0.9) 

D(m) 

Ps( W) 

Pr{ W) 

Ps( W) 

Pr( W) 

Ps( W) 

Pr( W) 

100 

0.00033 

0.0006 

0.000312 

0.000197 

0.000651 

0.00049 

200 

0.016 

0.0127 

0.0057 

0.0033 

0.0104 

0.0101 

300 

0.0877 

0.0595 

0.0289 

0.0223 

0.0527 

0.0525 

400 

0.3366 

0.1537 

0.0706 

0.0703 

0.1666 

0.0198 

480 

0.6225 

0.3601 

0.1466 

0.1455 

0.3455 

0.0342 

620 

- 

- 

0.4107 

0.402 

0.9627 

0.0373 

650 

- 

- 

0.4970 

0.483 

- 

- 



Fig. 7. Energy consumption per information bit versus source-destination distance when the relay is located in the middle 
over spatially correlated Nakagami—0.75 fading channels at target BER of 10 -2 for 4—QAM/QPSK constellation and different 
spatial correlation values. 


is assigned to the source and relay nodes in order to overcome performance losses incurred by 
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TABLE III 

Optimal Transmit Power Values for Source and Relay for relay located in the middle over correlated 
Nakagami— 1.25 Links with target BER of 10 -2 Using 4—QAM / QPSK Modulation. 


Corr. 

P = o 

p = 0.5 

p = 0.9 

D(m) 

Rs(W) 

Pr{ W) 

Rs(W) 

Pr{ W) 

Rs(W) 

Pr( W) 

100 

0.000372 

0.000197 

0.000424 

0.000233 

0.000514 

0.000296 

200 

0.0057 

0.0033 

0.0064 

0.0039 

0.0077 

0.0049 

300 

0.0225 

0.0223 

0.0249 

0.0249 

0.0290 

0.0290 

400 

0.0706 

0.0703 

0.0788 

0.0785 

0.0918 

0.0915 

500 

0.1728 

0.1712 

0.1929 

0.1912 

0.2249 

0.2228 

600 

0.3598 

0.3530 

0.4021 

0.3944 

0.4692 

0.4602 

630 

0.4381 

0.4281 

0.4897 

0.4785 

- 

- 

650 

0.4970 

0.4843 

- 

- 

- 

- 


the spatial correlation, within the given resource constraints. Concrete examples are shown in 
Table III and Table IV for some indicative transmission distances and the two target BERs; there, 

_ (J _ JJ 

the energy savings using CT scheme, defined as 1 — E /E for p = 0.9, p = 0.5 and p = 0 for 
the target BER of 1CT 2 at ds.u = 390m are 65%, 67% and 69.3%, respectively, whereas for a 
target BER of 10” 3 at ds,D = 250m the energy reduction is 69%, 72% and 73.8%, respectively. 
Interestingly, beyond a critical distance of 320m, even the highly spatially correlated CT mode 
at target BER of 10' 3 exhibits better energy-efficiency than the DT scheme with target BER of 
10 ~ 2 . 

In the same context, Fig. 7 illustrates the average total energy per information bit for both direct 
and cooperative transmissions for m = 0.75, which corresponds to severe fading conditions. The 
target BER is set to 1CT 2 under the given transmit power constraint while p = {0, 0.5, 0.9}. The 
transmit power is again allocated optimally to the source and relay with the latter positioned in 
the center of the network as shown in Table V. It is shown that DT outperforms CT only when 
ds,D < 170m. However, as the distance increases beyond this point, the corresponding overall 
benefits by CT are significant even under the worst spatial correlation scenario. Indicatively, 
at a transmission distance of 280m, which is the maximum distance that DT can operate with 
the available maximum transmission power, the energy gains by CT are 68%, 72% and 74%, 
respectively. In addition, it is shown that the advantage of the cooperation is more significant in 



















32 


TABLE IV 

Optimal Transmit Power Values for Source when the Relay is Located in the Middle over Correlated 
Nakagami—1.25 Links with Target BER of 10~ 3 Using 4-QAM / QPSK Constellations. 


Corr. 

P = o 

p = 0.5 

p = 0.9 

D(m) 

Ps{ W) 

Pr{ W) 

Ps{ W) 

Pr{ W) 

Ps{ W) 

Pr{ W) 

100 

0.00012 

0.0007 

0.00014 

0.0009 

0.00022 

0.00014 

200 

0.0179 

0.0116 

0.0218 

0.0147 

0.0331 

0.0233 

300 

0.0762 

0.0762 

0.0914 

0.0914 

0.1332 

0.1332 

410 

0.2662 

0.2651 

0.3196 

0.3183 

0.4660 

0.4641 

440 

0.3533 

0.3513 

0.4242 

0.4218 

- 

- 

480 

0.5009 

0.4969 

- 

- 

- 

- 


TABLE V 

Optimal Transmit Power Values for Source and Relay when the Relay is Located in the Middle over 
Spatially Correlated Nakagami—0.75 Links with target BER of 10~ 2 Using 4-QAM / QPSK Modulations. 


Corr. 

P = o 

p = 0.5 

p = 0.9 

D(m) 

Ps( W) 

Pr( W) 

Ps( W) 

Pr( W) 

Ps( W) 

Pr( W) 

100 

0.000754 

0.0004635 

0.0009112 

0.0005811 

0.00013 

0.0009 

200 

0.0119 

0.0076 

0.0143 

0.0093 

0.0386 

0.0097 

300 

0.0509 

0.0509 

0.0608 

0.0608 

0.0873 

0.0873 

400 

0.1612 

0.1606 

0.1926 

0.1918 

0.2765 

0.2753 

460 

0.2823 

0.2802 

0.3373 

0.3348 

0.4845 

0.4808 

500 

0.3945 

0.3905 

0.4714 

0.4666 

- 

- 

530 

0.4985 

0.4921 

- 

- 

- 

- 


severe fading conditions. 

Fig. 8 illustrates the average total energy consumption per information bit required for CT and 
DT as a function of S-D distance for fading parameters, m = {0.75,1.25,1.75, 2.25} for CT and 
m = 2.25 for DT. The target BER is set to 1CT 3 , the transmit power is allocated optimally, the 
relay is located in the middle and a zero spatial correlation is assumed. It is observed that the 
critical distances below which DT outperforms the CT in terms of energy efficiency are 250m, 
230m and 220m for m = 1.25, m = 1.75, and m = 2.25 , respectively. Moreover, the analysis 
indicates that DT with non-severe multipath fading condition (Nakagami—2.25) can operate only 
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x 10“ 6 



Fig. 8. Energy consumption per information bit versus source-destination distance when the relay is located in the middle over 
uncorrelated Nakagami —m fading conditions at target BER of 1CP 3 for 4—QAM / QPSK constellation. 


up to 360m before utilizing the maximum transmit power, while the CT extends significantly 
even for moderate fading conditions, except for the worst case scenario (m = 0.75). It is also 
shown that the gain from the cooperation is not uniform as the Nakagami parameter increases 
from m = 0.75 tom = 1.25, from m = 1.25 to m = 1.75 and then from m = 1.75 tom = 2.25. 

Finally, Fig. 9 depicts the cooperation gain, defined in (l4TT) . when the relay is located in the 
middle of the source and destination. The power allocation is again carried out by using the 
derived OPA scheme for target BER of 10~ 2 under the maximum transmission power constraint 
with p = {0, 0.5, 0.9} for m = 1.25. The transmission distance is limited to 390m since beyond 
this limit it is only the CT mode that can transmit until its maximum transmission distance, 
depending on the spatial correlation between S-D and R-D paths. When the cooperation gain 
is below unity, the DT is actually more energy efficient than CT. As already mentioned, the 
reason behind this is that when CG < 1, which corresponds to relatively small transmission 
distances, the actual transmit power constitutes only a small fraction of the total average power 
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Fig. 9. Cooperation gain versus source-destination distance when the relay is located in the middle over spatially correlated 
Nakagami—1.25 fading environment at target BER of 1CU 2 for 4—QAM / QPSK constellations. After source-destination distance 
of 390m, direct transmission runs out of power and cannot anymore reach the target BER. 


consumption. However, when CG > 1, the system benefits significantly from cooperation and 
in general CG increases proportionally as the transmission distance increases for all scenarios 
of spatial correlation between the S-D and R-D paths. Interestingly, the existence of such 
efficiency threshold distance also implies that a hybrid system, where cooperation is only sought 
and deployed beyond certain minimum distance, can provide the most comprehensive solution 
to the energy-efficiency optimization. The analysis and modeling results and tools provided 
in this article form directly the basis for further development of such schemes in different 
communication scenarios, which forms an important topic of future work. 

VIII. Conclusions 

This work was devoted to the end-to-end SER analysis as well as the energy efficiency analysis 
and optimization of both direct and regenerative cooperative transmissions over Nakagami —m 
fading conditions in the presence of spatial correlation. Novel closed-form expressions were 
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firstly derived for the symbol-error rate of both M —QAM and iff— PSK constellations which 
were subsequently employed in formulating the constrained energy analysis and optimization 
problems under destination bit-error-rate target and maximum transmit power constraints con¬ 
sidering both transmit energy as well as the energy consumed by the transceiver circuits. The 
corresponding results indicate that depending on the severity of multipath fading, spatial corre¬ 
lation between the source-destination and relay-destination paths and the location of the relay 
node, the direct transmission can be more energy efficient than cooperative transmission but only 
for rather short transmission distances and up to a certain threshold value. Beyond this value, the 
system, as expected, benefits substantially from relaying and the corresponding cooperation gain 
increases proportionally to the transmission distance. It is expected that the offered results can be 
useful in the design, dimensioning and deployment of low-cost and energy efficient cooperative 
communication systems in the future, especially towards the green communications era where 
the requirements and incentives towards energy consumption optimization are considered critical. 


Appendix 


A. Proof of Theorem 1 


The average SER for the direct transmission can be expressed as, 



(61) 


Evidently, a closed-form expression for (RTTt) is subject to analytic evaluation of the following 
integrals 



(62) 


and 



(63) 


By re-writing the indefinite form of the above class of integrals as 



( 64 ) 
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and setting u = cos 2 {0), one obtains 


Z(a, m) = — 


(i -y) r 


-d u. 


(65) 


2y/u(l — u + a) r 

The above integral can be expressed in closed-form in terms of the Appell hypergeometric 
function of the first kind, namely 


cos(6) (11 o 2 

M a , m) = , ^ m F 1 ( - m, m, cos (0), 


3 

2 ’ 


cos 2 (6) 


(1 + a) m ~ ^ \2 J 2 2’"~" 1 + a 

Equation (l66l) reduces to zero when 9 = tt/ 2. To this effect, it immediately follows that 


( 66 ) 


1 (a,m, 0,^ = 


PI 11 3 i 


(l + a) m_i V 2, 2 2 , “’1 + a, 

which with the aid of the properties of Ti(.) function can be equivalently expressed as 


(67) 


1 (a,m, 0, 

In the same context, 


V^T (m + I) 


. , . , oFi ( m, m + 1 ,- 

2(1 + a) m T(m + 1) V 2 1 + a 


( 68 ) 




PI 11 3 1 1 

1 2’ 2 ’ ’ 2’ ’ 1 + a 


1 A 1 3 1 1 

—=Ti -m, m, -, -, —-r 

y/2 v 2 2 2 2 2(1 +a) 

(69) 


which can be alternatively expressed as follows: 


I[a,m, 0, — ) = 


(m + 2 F 1 + yA) F 1 \ m,m, §, 2( i+ a) ) 


(70) 


2(1+ a) m T(m +1) ^2(1 + a) m 

By performing the necessary change of variables in (1681) and (T70l) and substituting in (|6T1) yields 

©. 


B. Proof of Lemma 1 

The J (a, b, m) integral can be re-written as 




sin 4m (0) 

(sin 4 (6>) + asin 2 (6 l ) + b) m 


d9. 


(71) 
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By setting u = cos 2 (6) it follows that 


J(a,b,m) = —- 


(l-«) 


2 m 


2 J y/uy/1 — u [(1 — u) 2 + a( 1 — u) + 6]’ 

By applying the binomial theorem in lT74l eq. (1.111)] one obtains 


-da. 


(72) 


2m-§ 

J(a,b,m) = - £ (2m ~ 2 
1=0 


‘I)' 


l— I 
U 2 


-du. (73) 


l J 2 J [1 — u(a + 2) + u 2 + a + b}' 

The above integral can be expressed in terms of the Appell function of the first kind yielding 


c lm—\ 


J(a,b,m ) - - y. 


1=0 

x(l 


2m 


i\ (-l)^M i+ 2 Fi ^ l + |, 


2n 


2n 


2+a—Va 2 —46 ’ 2+a+\/a 2 —46 


2m 


2 + a — \/a 2 — 46 


(1 + 2/) (a + b + (w — l) 2 — aw) f 
2 u 

2 + a — \/a 2 — 46 


(74) 


which upon performing the necessary counter-substitution and algebraic manipulations yields 


C. Proof of Lemma 2 

The JC(a,b,m,n) integral can be alternatively re-written as 


/C(a, 6, m, n ) = 


Sill 


2m+2n 


(*) 


(sin 2 (6 l ) + a) m (sin 2 (6 l ) + 6)' 
which by setting u = cos 2 (6 l ) can be expressed as 


-d 9 


(75) 


/C(a, 6, m, n) = 


(1-m) 


m+n—i 


dw. 


2(1 - M + a) m (l - w + 6) ? 

By applying the binomial theorem in lf74l eq. (1. Ill)], it immediately follows that 


(76) 


m+n— ^ 

/C(a, 6, m, n) = 

1=0 


V I 


(-i y 

2 


u 


i- i 


(1 — m + a) m (l -m+6) 5 


-da. 


(77) 


The above integral can be expressed in closed-form in terms of the Appell function of the first 
ki nd. As a result, by making the necessary counter-substitution and performing some long but 
basic algebraic manipulations, equation (fTQl) is deduced. 
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D. Proof of Lemma 3 

By setting u = sin 2 (6 l ), it follows that 


/ 


sin 2m (0)d0 


/ 



(78) 


The above integral can be expressed in closed-form in terms of the Gaussian hypergeometric 
function, namely, 



(79) 


Therefore, by performing the counter-substitution, equation (fl5l) is deduced. 

E. Proof of Convexity of the Optimization Problem 

Below we prove the existence of optimal powers, which are subsequently employed in min¬ 
imizing the overall energy consumption in the considered cooperative communication system. 
Based on (l39l) . the average total power consumption can be re-written as follows: 


Pt = (Ci + (1 + a)P s ) + (C 2 + (1 + a)P R ) (1 - SER s ,r) 


(80) 


where C\ = P CTr + 2P C R x and C 2 = Pct x + Pcr x • The symbol-error-rate SER, s r can be 
expressed in closed-form with the aid of Theorem 1. This expression is a function of 1 + 
P s Cl S p Uqam/ Lf s R N 0 rris.nl therefore, for proving the existence of the optimum values, it is 
sufficient to show that () 2 P^'/ if 2 If > 0 and <f 2 P ( ' r j(f 2 P 2 > 0. To this end, it is straightforward 
to show that <9 2 P^/<9 2 P| = 0. Likewise, based on the optimal condition in (l60l) and taking the 
second-order partial derivative w.r.t the variables Ps and P R , one obtains 


8 2 Pt d 2 P° T 


d 2 p2 - Q2p s p R 


(81) 


where 


d 2 Pf 

dP s P R 



(1 + a i)N 0 P Ls , R 


(82) 


and 
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K a = 


4C' / P s n s , R ()QAM 40 2 / PgClg,R 9 qam 

— 1 W~5 - ' ms ’ R ' o - 1 - ,™s,R, 

n \N 0 P LsR m s ,R 2 ) tt \N 0 P LsR m s , R 


(83) 


where it is recalled that C = 1 — 1 /\[M for the case of M—QAM whereas A' 4 denotes the SER 
representation with values in the range 0 < SER < 1 and with all other constants being positive. 
Based on this, the second-order partial derivatives w.r.t Pg and P R are always greater than or 
equal to zero, which implies that d 2 P T /d 2 Pg > 0. Given the general second order conditions 
in IfTTTl . it immediately follows that (l80l) is convex w.r.t to Pg and P R and possesses a unique 
minimum value. 
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